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Abstract 

We assume the cosmological dark sector to consist of pressureless matter and holographic dark 
energy with a cutoff length proportional to the Ricci scale. The requirement of separate energy- 
momentum conservation of the components is shown to establish a relation between the matter 
fraction and the (necessarily time-dependent) equation-of-state parameter of the dark energy. Sin- 
gling out intrinsically adiabatic pressure perturbations of the dark-energy component, the matter 
perturbations are found as linear combinations of the total energy-density perturbations of the 
cosmic medium and the relative (non-adiabatic) perturbations of the components. The resulting 
background dynamics is consistent with observations from supernovae of type la, baryonic acoustic 
oscillations and the differential age of old objects. For a physically interesting parameter space the 
perturbation quantities show a stable behavior. On smaller scales the model predicts oscillations 
in the dark sector which potentially should be testable by weak gravitational lensing and by their 
impact on the anisotropy spectrum of the cosmic microwave background. 
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I. INTRODUCTION 



Despite of the many efforts to consistently explain the results of the observations of 
supernovae of type la (SNIa) in the physical nature of the dark sector of the Universe, 
assumed to consist of dark matter (DM) and dark energy (DE), remains largely mysterious. 
The favored cosmological model is the A-cold-dark-matter (ACDM) model which also serves 
as a reference for alternative approaches to the DE problem. Grosso modo, the ACDM 
model does well in fitting most observational data (see, e.g., the recent WMAP 9 results 
Nevertheless, there is an ongoing interest in alternative models within General Relativity 
(GR) itself and beyond it. Although no serious contender seems to be around at the present 
time, these efforts continue to make sense not only because of the notorious cosmological 
constant and coincidence problems (see, e.g. but also to test as many potential deviations 
from the "standard" description as possible in order to constrain additional parameter sets 
which are usually introduced in alternative approaches in order to quantify these deviations. 
Typically, these approaches "dynamize" the cosmological constant in terms of scalar fields 
or fluids with a generally time-dependent equation of state (EoS). Constraining a potential 
time dependence by the data, e.g., is then crucial for a comparison with the ACDM model. 
Holographic models represent a specific class of dynamic approaches to the DE problem. 
These models are characterized by a relation between an ultraviolet cutoff and an infrared 
cutoff {4-6]. Such relation guarantees that the energy in a given volume does not exceed the 
energy of a black hole of the same size. The infrared cutoff has to be a cosmological length 
scale. For the most obvious choice, the Hubble radius, only models in which dark matter and 
dark energy are interactingwith each other also non-gravitationally, give rise to a suitable 
dynamics p, Q . Following [5], there has been a considerable number of investigations based 



on the future event horizon as cutoff scale 



9j. However, all models with a cutoff at the 



future event horizon suffer from the serious drawback that they cannot describe a transition 
from decelerated to accelerated expansion. A future event horizon does not exist during the 
period of decelerated expansion. A further option that has received attention more recently 
and which will be the subject of the present paper is a model based on a cutoff length 
proportional to the Ricci scale. The role of a distance proportional to the Ricci scale as 
a causal connection scale f or p erturbations was noticed in [10|. As a cutoff length in DE 
models it was first used in |ll| . Further investigations along this line have been carried out 
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in 12[. In 13Ml5| observational constraints were obtained on the basis of which Ricci DE 



was compared with the ACDM model. The dynamics of perturbations was considered in 



16 



181 ] . A number of studies performed in 19M24| include interactions between DM and 



25]. 



Ricci DE. A relation to quantum field theory has been claimed in 

In the present paper we reconsider the dynamics of a two-component system of pres- 
sureless DM and Ricci-type DE both in the homogeneous and isotropic background and on 
the perturbative level. While most dynamic DE scenarios start with an assumption for the 
EoS parameter for the DE, the starting point of holographic models is an expression for the 
DE energy density from which the EoS is then derived. Moreover, as was pointed out in 



261 ]. the mere definition of the holographic DE density generally implies an interaction with 
the DM component. Requiring this interaction to vanish imposes an additional condition 
on the dynamics. Non-interacting Ricci-type DE in particular is characterized by a simple 
relation between the matter fraction and the (time-dependent) EoS parameter. We shall 
confront the resulting background dynamics with recent SNIa data, results from baryonic 
acoustic oscillations (BAO) and from the history H(z) of the Hubble parameter. Based on 
a gauge-invariant analysis, the matter perturbations are found as a combination of the total 
and the relative energy- density perturbations. 

The structure of the paper is as follows. In Section [TT] we recall basic relations for holo- 
graphic models of DE. The resulting homogeneous and isotropic background dynamics is 
confronted with observational data in Section IIHI Section IIVI provides us with the general 
two-component dynamics of the cosmic medium. The first-order perturbation theory of the 
model is presented in sections |V] |VT] and IVHI On this basis, the final set of coupled equa- 
tions for the non-adiabatic perturbation dynamics is found in section IVIIIl A summary of 
the paper is given in the final section IIXI 

II. BACKGROUND DYNAMICS FOR RICCI DARK ENERGY 

We start by recalling basic features of holographic DE models in a homogeneous and 



isotropic background [26|]. The cosmic medium is assumed to be describable by pressureless 
DM with energy density p m and a holographic DE component with energy density pu- In 
the spatially flat case Friedmann's equation is 

3H 2 = 87rG(p m + p H ) . (1) 



In general, both components are not necessarily conserved separately but obey the balance 
equations 

p m + 3Hp m = Q , p H + 3H(1 + u)p H = -Q , (2) 

such that the total energy p = p m + Ph is conserved. Here, uj = = is the equation- 
of-state (EoS) parameter of the DE and pn is the pressure associated with the holographic 
component. The acceleration equation can be written 



H = — H 2 1 + 



uj 



1 + r 



d\nH 3 / uj{a) 
d\na 2 V 1 + r(a) 



(3) 



where r = — is the ratio of the energy densities. The total effective EoS of the cosmic 
medium is 



According to the balance equations (j 

f = 3Hr (1 + r 



o 1 + r 

the ratio r changes as 
uj Q 



Following (4, 5], we write the holographic energy density as 



1 + r + 3Hp, 



Ph 



3 c 2 M 2 p 
L 2 



(4) 



(5) 



(6) 



The quantity L is the infrared (IR) cutoff scale and M p 



8n G is the reduced Planck 



mass. The numerical constant c 2 determines the degree of saturation of the condition 

L 3 p H <M 2 Pl L, 



(7) 



which is crucial for any holographic DE model. It states that the energy in a box of size L 
should not exceed the energy of a black hole of the same size ^]. 

Differentiation of the expression (J5]) and use of the energy balances fl2]) yields 



Ph L 



(8) 



In general, there is no reason for Q to vanish. Assuming Q = provides us with a specific 
relationship between u and the ratio of the rates ^ and H. Any non- vanishing Q will modify 
this relationship. 

With Q from (jSJ), the general dynamics (E]) of the energy density ratio r becomes 

2 L 



-3H 1 



uj 



1 + r 3 HL 



(9) 
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The case without interaction is characterized by (cf. eq. (J5])) 



Q = Q f = r [ 2j -3H \ = 3Hru 



(10) 



with a generally time-dependent u. Different choices of the cutoff scale L give rise to different 
expressions for the total effective EoS parameter in (jlj) and to different relations between 
u and r. Our interest in the present paper will be the Ricci-scale cutoff. The role of a 
distance proportional to the Ricci scale as a causal connection scale for perturbations was 



noticed in 



lOj. In 11] it was used for the first time as a DE cutoff scale. The Ricci scalar 



is R = 6 ( 2H 2 + H ) . For the corresponding cutoff-scale one has L 2 = Q/R, i.e., 

R 



p H = 3c 2 M 2 



where a = -^q. Upon using ([3]) we obtain 

a 



6 



a[2H z + H 



(11) 



PH 



— H 



UJ 



(12) 



2 V 1 + r, 

for the holographic DE density. Notice that the (not yet known) EoS parameter explicitly 
enters pu- Use of Friedmann's equation provides us with 



1 = - (1 + r - 3w) 



1 / 2 

u= - (1 + r) ~ 

3 V ; 3c 2 



(13) 



Obviously, a constant value of u necessarily implies a constant r and vice versa. The time 
derivatives of u and r are related by f = 3u. The second relation in (JT3"|) can be used to 
express c 2 in terms of the present values (subindex 0) of u and r: 



— = 1 + r - 3u>, 



o • 



r = r + 3(u- UJ ) . 



(14) 



The parameter c is related both to r and ojq. In a next step we differentiate pu in ( TTTj) 
which yields 

p H = a(AHH + h) . (15) 



With the help of (|3]) and the definition (TTT1) with we derive 26] 



p H + 3i/(l + u;)pH = -Q , 



where 



3iJ 



1 + r 



Ph ■ 



(16) 



(17) 



Relation (HM with (ITTj) . which implies that in the general case one has Q ^ 0, i.e., both 
dark components do interact with each other also non-gravitationally, is a direct con- 
sequence of the ansatz ( II ip . The DE balance in Eq. ([2]) may then also be written as 
Ph + 3H (1 + oj e ff) ph = with an effective EoS parameter 

The present ratio r is related to the present matter fraction Q m0 of the Universe by r = 

1 — f^mO ' 

According to relation (fl71) . a constant EoS parameter to is compatible with Q = only 
for co = 0, i.e., if behaves as dust. If we admit co ^ 0, however, there exists a non-trivial 
case Q = 0: 

^ w din a; 

Q = ra = - r = — . 19 

ii a ma 

It is this configuration that we shall investigate in the present paper. Equation (fl9|) together 
with the second relation of f TH|) is a differential equation for to which has the solution 

U = U \ a-(r -^-3oo ' (20) 
where we have normalized the present value of the cosmic scale factor a to ao = 1. There is 
no freedom left to choose the equation of state. It is fixed by the choice of pu together with 
the requirement Q = 0. Notice that this is different from the more familiar procedure to deal 
with (non-holographic) DE, where one starts with an assumption for the EoS parameter and 
afterwards finds an expression for the DE density by integrating the corresponding balance 
equation. Here, the starting point is the energy density and the EoS parameter has to be 
derived. 

Knowing the EoS parameter f )20|) . it follows from (TH|) that 

r - rn r ° ~ 3uj ° fan 

°r - 3^0-3-0) ' 

i.e., r(a) is fixed as well. At high redshifts we have 

to ->■ , r -)■ r - 3co , (a < 1) . (22) 



The property t 
pointed out in 



rat non-interacting Ricci-DE behaves as dust at high redshift was already 



11] . The values in the far- future limit are 



co -> co - ^r , r -> , (a > 1) . (23) 
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The limits in ( 122]) imply that this model naturally reproduces an early matter-dominated 
era. For ro ~ | and Uq ~ —1, the ratio r approaches r « y for a < 1. This value is only 
roughly ten times larger than the present value ro. For the ACDM model the corresponding 
difference is about nine orders of magnitude. In this sense, the coincidence problem is 
considerably alleviated for the the present model. On the other hand, in the opposite limit 
a ^> 1 the ratio r approaches zero as for the ACDM model. Apparently, the far-future EoS 
can be of the phantom type for uq — |ro < — 1. 
The Hubble rate of our model turns out to be 



^_ = a _ 3 /2 Su a^- 3 ^ -r [l + r -3uJo 
H V 3uj - r [1 + r - 3w ] 



(24) 



For a < 1 we recover the Einstein - de Sitter behavior H oc a 3 / 2 . The total effective EoS is 

P = g = r - 3u 

p 1 + r W ° r a-( r °- 3w °) [1 + r - 3u ] - 3u ' 

For the adiabatic sound speed of the DE component we find 



Ph V 3 1 + UJ 

and the corresponding quantity of the total cosmic medium is 

P = P Ji = J-^ . (27) 
P P 1 + T7^ 

With the solutions ( 120]) and (12ip all these quantities are explicitly known, i.e., the background 
dynamics is completely solved analytically. In the following section we perform an actualized 
confrontation of the background dynamics with recent observational data. 



III. OBSERVATIONAL ANALYSIS 



Our observational analysis of the background dynamics uses the following three tests: the 
differential age of old objects based on the H(z) dependence as well as the data from SNIa 
and from BAO. A fourth test could potentially be added: the position of the first peak of the 
anisotropy spectrum of the cosmic microwave background radiation (CMB). However, the 
CMB test implies integration of the background equations until z ~ 1.000 which requires 
the introduction of the radiative component. But the inclusion of such radiative component 
considerably changes the structure of the equations and no analytic expression for H(z) is 
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available. Hence, we shall limit ourselves to the mentioned three tests for which a reliable 
estimation is possible. 

Based on the evaluation of the age of old galaxies that have evolved passively, there are 



13 observational data available for the differential age 
data has been considered 32| . The basic relation is 



27H31|. Recently, a new set of 21 



ldz 

H(z) = — 

K J zdt 



(28) 



The value of the Hubble parameter today can be added to these data, leading to 14 or 22 
observational points, depending on the sample used. 

The SNIa test is based on the distance modulus /i which is related to the luminosity 
distance D L by 



M = 51og 10 D L 



D, 



dz' 



(29) 



H Jo ^/Hjz 7 ) ' 

In this expression we have restored the velocity of light c. The quantities m and M denote 
the apparent and the absolute magnitudes, respectively. 

Two decisions have to be taken for this test. The first one concerns the choice of the 
sample. There are many different SNIa data sets, obtained with different techniques. In 
some cases, these different samples may give very different results. The second point is the 
existence of two different calibration methods. One of them uses cosmological relations and 
takes into account SNIa with hi gh z (Salt 2), the other one, using astrophysical methods, is 
suitable for small z (MLCS2k2) 33]. In some cases, the application of different calibrations 
can lead to different results also. All this makes the SNIa analysis very delicate. Here, we 
use the Union 2 sample [34j, calibrated by the Salt 2 method. 

Baryonic acoustic oscillations have their origin in oscillations in the photon-baryon plasma 
at the moment of the decoupling at about z = 1.090. They can be characterized by the 
distance scale I35I , 

a .,^.,[1 r^r (30) 



JH(z b )]^[z b J H(z) 
We shall use the WiggleZ-data [36] A = 0.474 ± 0.034, 0.442 ± 0.020 and 0.424 ± 0.021 for 

the redshifts z b = 0.44, 0.60 and 0.73, respectively. 

Generally, the key quantity of a statistical analysis is the y 2 parameter 



a* 



(31) 



FIG. 1: One dimensional PDFs for ojq (upper panels) and ro (lower panels) using, from left to 
right, H(z), SNIa, BAO and the combination of the three tests. 

where f{x 3 )\ is the theoretical evaluation of a given observable, depending on x- 7 free param- 
eters, f{xi)° is the corresponding observational value with an error bar <jj and n is the total 
number of observational data for the given test. In terms of the \ 2 parameter one defines 
the probability distribution function (PDF) by 

P(x j ) = Ae~ x2/ \ (32) 

where A is a normalization constant. The estimations for one or for two given parameters 
are obtained by integrating over the remaining ones. For a combination of all tests we use 
the total x 2 -value Xt; 

Xt = Xh( z ) + X 2 sn + X 2 bao- ( 33 ) 

Assuming a spatially flat universe, the three free parameters of the model are the density- 
ratio parameter r = , , the EoS u and the reduced Hubble parameter h, defined by 
H = 100/ikm/s/Mpc. 

In figure [T] we display the one-dimensional PDFs for each of the tests and for their 
combination. The results for the density parameter are different for each test. The 
combination of all tests leads to a value of r ~ 0.4, corresponding to Vt m o ~ 0.29, roughly 
in agreement with the ACDM model. For the equation of state parameter we obtain 
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FIG. 2: Two-dimensional PDF for uq and ro resulting from a combination of the three tests. 



Uo ~ —1, consistent with the ACDM model as well. According to the first relation of 
( JT41) . the parameter c 2 turns out to be c 2 ~ 0.46. This value is coincides with the result 
in The two-dimensional PDFs at la (68% of confidence level), 2a (95% of confidence 
level) and 3a (99% of confidence level) are shown in figure [2] The estimation for uq, 



based on a combination of the three tests at 2a, is uo = —0.987! 



-0.083 
-0.100' 



while for ro we find 



r = 0.406! 



s-0.073 
-0.061- 



IV. GENERAL TWO-COMPONENT DYNAMICS 

To study the dynamics of inhomogeneities on the homogeneous and isotropic background 
of the previous sections we first consider the general description of the two-component sys- 
tem. It is characterized by a total energy-momentum tensor 

T ik = puiu k + ph ik , T^ k = 0, (34) 

where hik = Qik + UiUk and gikU l u k = — 1. The quantity u % denotes the total four- velocity 
of the cosmic substratum. Latin indices run from to 3. The total splits into a matter 
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component (subindex m) and a (holographic) DE component (subindex H), 

T ik = T ik + T ik^ ^ 

with (A = m, H) 

T k = PAu\u k A + p A K k , hf = g lk + u\u k A . (36) 
For separately conserved fluids we have 

T k ]k = 0, T k . k = . (37) 

Then, the separate energy conservation equations are 

u iniT m = p mya U m + QrnPrn = (38) 

and 

- u Hi Tf[ . k = p H ^u a H + &h (ph + Ph) = . (39) 

In general, each component has its own four- velocity, with gikU l A u k A = — 1. The quantities 
Qa are defined as = u A . a . For the homogeneous and isotropic background we assume 
u m = u h = M<1 - Likewise, the momentum conservations are written as 

^mi^m :k = PmU m = (40) 



and 



where ii A = u A . b u A . 



h a m n ± = (Ph + Ph) u a H + PhAh = 0, (41) 



Denoting first-order perturbations about the homogeneous and isotropic background by 
a hat symbol, the perturbed time components of the four-velocities are 

u = u° = u° m = u° H = \goo ■ (42) 



According to the perfect-fluid structure of both the total energy-momentum tensor 
and the energy- momentum tensors of the components in ( |36i) . and with = u% = u a in 
the background, we have first-order energy-density perturbations p = p m + Ph, pressure 
perturbations p = p m + Ph = Ph and 

= T^ a + f Ha => (p + p)u a = p m u ma + (p H + p H ) u Ha . (43) 
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In linear order the spatial components of the 4— accelerations are 

U^ = U^~ -g00,fM, U(rn)fA, = U( m )^ - -000,/^ W(H)^ = W(ff> ~ ^900,fi ■ (44) 

For the first-order pressure gradient terms we find (recall pu = p) 

P,ch c Hfl = p,» + pu H », P,ch% = p ijU + pu» . (45) 
From the matter-momentum conservation (1401) it follows that 



u* = =4> U( m j M - ~<?oo lAt = =^> M (m)At = ^oo,/. • (46) 

According to ( 143]) the differences between w Q and the corresponding quantities of the com- 
ponents are 

A , _ t — P m (" _ t \ - _ r — ^ H +P g / - _ - \ ci7^ 

P +P P +P 

Restricting ourselves to scalar perturbations, the perturbed line element can be written 
ds 2 = - (1 + 20) dt 2 + 2a 2 F a dtdx a + a 2 [(1 - 2f) 8 a p + 2E ta p] dx a dx? . (48) 
Furthermore, we define the (three-) scalar quantities v, v m and Vh by 

a 2 w M + a 2 F M = u lt = v ttl (49) 

and 

a 2 ^ + a 2 F^ = u mtl = v m ^ , a 2 u^ + a 2 F^ = u Hfl = v H ,^ , (50) 

respectively. 

From the definitions of O, G m and 0# it follows that 

B = (Av + Ax) - 3ip - 3H4> , (51) 
where A is the three-dimensional Laplacian, 

x = a 2 (E-F) (52) 

and 

6 m = \ (Av m + A X ) - 3^ - 3H<P , 6 H = ^ + A X) - 3ip - 3#0 . (53) 
cr er 
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V. CONSERVATION EQUATIONS IN LINEAR ORDER 



At first order, the energy balances (1381) and (|39l) are 

Pm + PmU° + Q m Pm + @P m = (54) 

and 

'pn + Phu° + &h (ph + Ph) + ® (ph +Ph)=0 , (55) 
respectively. The total first-order energy conservation takes the form 

'p + pu° + 6 (p + p) + 6 (p + p) = . (56) 

Comparing (J3i|) . fl55|) and fl56l) one finds 

e = ^e m + (57) 
p+p p+p 



with p = ph- The separate momentum conservation equations are given by poj) and (1411) . 
Additionally, we have the total momentum conservation 

(pm + pH+p H )u a +PH,ih ai = . (58) 

For the total energy-density perturbations S = £ equation ( 1561) yields 

j_e^5-^ + e + = , (59) 

p p \ pj p 

where we have used u° = —<f). 

The momentum balance (1551) for the cosmic medium as a whole together with the first 
relation of (1441) and the second relation of (1451) provides us with 

(p + p) [v + 4>] + P c = , (60) 

where p c = p + pv. 

In terms of the fractional perturbation 5 m = — , the matter energy conservation (|54p can 
be written as 

L-— + e m = o. (6i) 

Pm 

The matter momentum balance (140]) together with the second relation of ( 14"4"|) results in 

i> m + <P = . (62) 
13 



With 5h = — the energy conservation flo^]) for the DE component is 

5 H + p JL SH -PJL^ + e H (i + P») + G(8 H + i»)=o. (63) 

PH PH V PHj V PHj 

In the following it will be convenient to introduce the quantity 

D H = -^— = -^. (64) 
Ph+Ph 1 + u 



In terms of Dh eq. ( 16311 then takes the form 

Dh-Q^Dh + O^ + Qh + Q Ph =0. (65) 
Ph + Ph 

The dark-energy momentum balance ( HTj) together with the third relation of ( |4~4|) and the 
first relation of f )45|) result in 

(Ph+Ph)[^h + 0]+P c " =0 =► ^ + = ^^ = -P c h, (66) 

Ph + Ph 

where p c a = p H + p#t># = p + pt>//. 

Our final goal in this paper is to calculate the matter-density perturbations. To this 
purpose we shall solve the coupled system of total energy perturbations 5 and relative 
energy-density perturbations 5 m — Dh- In the following section we start by establishing the 
equation for the total energy-density perturbations. 

VI. TOTAL ENERGY-DENSITY PERTURBATIONS 

We consider equation ( 1591) and introduce therein the gauge-invariant quantities 

5 C = 5 + to , 6 C = 6 + Ov . (67) 

Then, eq. ( 159]) is rewritten as 

5 C - B-5 C - £ (v + <p) + -p c + ( 1 + -) B c = . (68) 
P P P V PJ 

Combination of the energy conservation (1681) and the momentum conservation (|60|) yields 

5 C - Q-5 C + ( 1 + ^ j 9 C = . (69) 



P V P, 

The perturbation B has to be determined from the perturbed Raychaudhuri equation for 
6. Neglecting shear and vorticity, the Raychaudhuri equation is 

<d + \o 2 -u% + 4nG(p + 3p) = . (70) 
14 



In terms of the gauge- invariant variables one finds, at linear order, 



e c + -ee c + a-kG P 5 c - < = o 



(71) 



In a next step we have to differentiate (1691) and to insert (ITT]) into the resulting expression. 
The remaining C terms can be eliminated by f )69|) again. Using also 

1 / Ap + pAv \ 1 Ap c 



u: 



P + P 



a 2 p + p 



(72) 



the equation for 5 C becomes 



5 C + 



2-6 P - + 3 P - 
P P. 



H5 C 



p 9p 2 



V 



P 2p2 



p. 



2 re 



1 Ap c 
a 2 p 







(73) 



Changing to a as independent variable (5' = and transforming into the k space, we 
arrive at 



5 C " + 



3 _ 15 p 3 P 
2 2 p p 



8«_ 
a 



3 p 9 p 2 p~\ 5 C k 2 p c 
l_ 12— - 9— | - 

2 pip 2 pa 2 a 2 H 2 pa 2 



. 



(74) 



Equation (|74j) governs the behavior of the total energy- density perturbations. As we shall 
see, via the p c term the perturbations S c are coupled to the relative perturbations 5 m — Dh- 



VII. COMBINING THE SEPARATE CONSERVATION EQUATIONS 



Now we combine the separate energy conservation equations flBT]) and fl65l) of the compo- 
nents and define S m u = S m — -D//. Then 



S m H — © ( -Pff _ T-^il ) + @ff — @m 



(75) 



where Ph = P ^" PH ■ Combining the momentum balances ( 1621) and ( 1661) results in 



(76) 



Because of the structure of the first-order expressions in ( l53|) one has 



©H - @m = —A (Vjj ~ V % 

a 



(77) 



Eq. fl75|) then becomes 



= O ( P H - ^-D H ) + ^jA (v H - v r , 
Ph J a 



(78) 
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Differentiation of ( 178]) yields 



S, 



mH 



PH 



- 2h\a (v h - v m ) + \a (v h -v m )' . 
a a 



Using here ( 178]) again and also (175]) results in 



PH 



Ph J a 2 



(79) 



(80) 



The difference Ph — f^^H describes the deviation of the DE pressure perturbations from 
being adiabatic. It is zero for purely adiabatic DE perturbations. We discuss this issue in 
more detail in the following section. 



VIII. NON-ADIABATICITY AND FINAL SET OF EQUATIONS 



Generally, the deviation from adiabaticity in a two-component system with components 
m (pressureless) and H is 

P VP _ pc _ P_ D c = Ph +Ph ( l>n i>n fit 



p + p p p + P 



P 



P + P \Ph+Ph PhPh+Ph 



Pm (ph +Ph)Ph 



(p+pY 



Ph 



Ph 



Ph+Ph 



Pm 
Pm 



Let us consider the combination p H — ^-pn- With p H = upn one has 

PH 



Ph . . / Ph 
Ph - —Ph = UJPH + [u - — 
Ph V Ph 



Ph 



where u is given by the solution ( [20]) and the adiabatic DE sound speed by 
of (1261) the combination ([82]) then results in 



Ph . 
Ph - —Ph 
Ph 



Ph 



(81) 

(82) 
Because 

(83) 



which is a gauge-invariant expression. In general, now an assumption for the perturbed EoS 
parameter to is necessary to proceed. We shall consider here the adiabatic case 



Ph . 
Ph = —Ph 
Ph 



to 



-~D H 

3 



(84) 



This assumption of an adiabatic DE component allows us to relate the otherwise undeter- 
mined perturbation to of the EoS parameter to the DE energy perturbation D^. Under this 
circumstances equation (IHUj) reduces to 

AP CH 



SmH + 2HS m H 
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(85) 



We emphasize that the total perturbation dynamics remains non-adiabatic due to the two- 
component nature of the medium. With an adiabatic DE component, the general relation 
(1ST]) simplifies to 



Pm (pH +Ph) PH 

p + P PP+P (p + pf PH 



s, 



mH 



or 



~c P ~c Pm (PH + PH) PH q 

P — ~P 7 : x —JraH 

P KP + P) PH 



P c = ? [p c - p m S mH ] 



(86) 



(87) 



Through f lH7|) the dynamics of the total energy- density perturbations, described by Eq. (1741) . 
is coupled to S m u- Explicitly, 



3 15 p p 

+ 3- 

2 2 p p 



6^ 
a 



- + 12^ - - 9^ - Jf_t 
2 P 2p 2 p a 2 H 2 p 



k 2 pp m S mH 
a 2 a 2 H 2 p p a 2 



At high redshift, for a< 1 equation (188]) approaches the Einstein - de Sitter limit 

(a < 1) . 

To obtain an expression for the term on the right-hand side of Eq. (1551) we write 



+ _ JL = o 

2a 2a 2 



p c H = p -f p-y^ = p + p<y -f p (y H — UJ 



(88) 



(89) 



(90) 



or 

p c H = p c -f p (^ _ w ) . 
Now, for the difference v — v h we have 



P + P 



(V m ~ V H ) 



From relation (J7SJ it follows that (in k space) 



V m — Vh — JpSmH 



Then p c n is written 



-c ~c • P m n 

P H =P -P^—T^SmH 

p + p k A 

Here we introduce p c from f )B7|) to obtain 

P 



p K 



p = - 

p 



p5 c — p m \ S m n — 3H—S m H 



(91) 



(92) 



(93) 



(94) 



(95) 
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By relation (195|) the dynamics of S m u in Eq. ( 150"]) is coupled to the dynamics of 8 C . Explicitly, 



D mH ' 



r p 



3 --3- 
2 1 + up 



3 p 
2 P. 



ql 



mH 



+ 



r p k 2 S. 



mH 



1 + up a 2 i/ 2 a 2 



1 + 



p k 2 5 C 
1 + uj J p a 2 H 2 a 2 



For a < 1 this equation approaches 



S" 



n,u. + 2a^ mH 







fa < r 



(96) 



(97) 



There is one decaying solution and a solution S m n = const. At high redshift, according 
to (1221) . the DE behaves as dust as well and we may use approximately adiabatic initial 
conditions S m H ~ for the coupled system (1551) and fl9"61) . 

Our interest is the matter-energy perturbation. To this purpose we decompose the total 
energy-density perturbation 5 C according to 



j-c Pm j-c . PH j-c 

° = °m + "H 

P P 



(98) 



Combination with S m H = 8 n 



4r~ leads to 



8° 



1 + 



1+r 



+ — JmH 

1 + r 



(99) 



which describes the matter-energy perturbations as a combination of 5 C and S m H- To obtain 
its dynamics one has to solve the coupled system of equations ( l8~8"j) and (|96|) . 

The result for the matter perturbations of our model can be compared with the behavior of 
matter perturbations in the ACDM model. The latter can be obtained as a limiting case from 
Eq. ( 1551) for the total energy- density perturbations. Using there p = Pa = —p\ = constant 
and 

1 



V 
P 



PA 

Pa + Pm 



"3 


15 p 


5 C ' 


2 ~ 


Tp_ 


a 



P = Pm + Pa , 
Eq. (1551) reduces to 
5 C " + 

For the perturbations we have 

p c = p c m =► & '- — 

P PA 

Then, the dynamics of 8^ is governed by 

3 2 + r 51 3 r 
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3 +12 P_9F 



p 2 p 2 \ a 
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(ACDM) 
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+ 
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2 1 + r a 



2 1 + r a 2 
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(ACDM) 



(100) 



(101) 



(102) 



(103) 



While for r 1 at high redshift the behavior of 8^ for our Ricci-type DE model is indistin- 
guishable from that of the ACDM model, the future evolution is different. The "growing" 
mode of the latter approaches 5^ = constant for r — > 0, independently of the scale, whereas 
for Ricci-DE the behavior of the corresponding quantity depends on the scale and decays 
and/or oscillates. It remains open, whether or not these (unwanted) oscillations will dis- 
appear if intrinsic non-adiabatic DE perturbations are included. This would amount to a 
modification of relation f )84|) . 

A final comment concerning the interpretation of the perturbation variable is in order 
here. The matter density perturbation in relation f[9"9~j) is defined with respect to the 
total comoving gauge. To obtain the matter density perturbation, comoving with the matter 
velocity, = 5 m + —v m , we have to consider 



Pn 



Since 



5 c m = 5 m + ^v = 6 m + ^v m +^(v-v m ) = 5 c - + — (v-v m ) . (104) 



Ph+Ph ( v a 2 - 

v-v m = ■ {v H -v m ), v H -v m = -—S mH , (105) 

p + p k 2 



the quantity of interest is 
5 



1+ 1- 



1 + r 



1 + t 1 - k 2 



u2 ' T \s' mH - (106) 



+r 

Obviously, b c m and 5^ differ by the last term in (11061) . Because of the factor ^M^- one expects 
that on scales smaller than the Hubble scale the differences between b c m and are small. 

The behavior of the quantities S m n, S c , b c m and 5^ together with 5^ of the ACDM model 
is visualized in Figures [3] and [5] for uq = —0.9 with k = 0.1 and k = 0.01, respectively. The 
parameters in Fig. H] differ from those of Fig. [3] by a higher energy- density ratio. This indi- 
cates a weak dependence on tq. All these perturbations show a stable behavior. Instabilities 
(not shown here) occur, however, for very small values of and for ujq < —1. 

All figures show a certain compensation effect between ripple-type features in S m u and 
5 C such that the matter perturbation curves b c m and look rather smooth. The figures 
also confirm that differences between b c m and 8^ are indeed small on the chosen scales. 

In figures [3] and H] there appear oscillations of all the perturbation quantities very close 
to the present time. They seem to be similar to those known from (generalized) Chaplygin- 
gas models. Potentially, such oscillations should be testable by weak gravitational lensing 
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0.0 0.2 0.4 0.6 0.8 1.0 

a 

FIG. 3: Evolution of the perturbation quantities S m jj (dotted line), 5 C (dash-dotted line), 5^ 
(dashed line) and 6%? (thin solid line) for loq = —0.9, ro = 0.4 and k = 0.1. For comparison the 
ACDM result (thick solid line) is also included. 



and by their possible influence on the CMB. But they may also jeopardize the model (cf. 
371 . l38j for a discussion of the matter power spectrum in the context of Chaplygin-gas 
models). On the other hand, we have chosen here the simplest possible case, assuming the 
intrinsic perturbations of the DE component to be adiabatic. It might well be that this 
is an oversimplification and a more general model that includes intrinsic non-adiabaticities 
might be more realistic. Moreover, the observed matter-spectrum refers to the distribution 
of baryons which so far are not included in our calculation. 



IX. CONCLUSIONS 



By its definition, non-interacting Ricci-type DE is characterized by a necessarily time- 
dependent EoS parameter. This makes it an observationally testable alternative to the 
ACDM model. Moreover, it establishes a a relationship between this EoS parameter and the 
matter content of the Universe. Ricci-DE behaves almost as dust at high redshift. The ratio 
of the energy densities of DM and DE varies considerably less than for the ACDM model. 
Since the time of radiation decoupling it has changed by about one order of magnitude 
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0.0 0.2 0.4 0.6 0.8 1.0 

a 

FIG. 4: The same quantities as in Fig. [3] for ojq = —0.9, tq = 0.8 and k = 0.1. 




0.0 0.2 0.4 0.6 0.8 1.0 

a 



FIG. 5: The same quantities as in Fig. [3] for loq = —0.9, ro = 0.4 and k = 0.01. 

compared with roughly nine orders of magnitude for the ACDM model. This amounts to a 
remarkable alleviation of the coincidence problem. Our statistical analysis, based on recent 
observational data from SNIa, BAO and H(z), results in a preferred value of c 2 ~ 0.46 
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for the Ricci-DE parameter which confirms earlier studies in the literature [ill]. Within a 
gauge-invariant analysis we calculated the matter perturbations as a combination of the total 
energy perturbations of the cosmic medium and the relative perturbations of the components 
and clarified the relation to the matter perturbations of the ACDM model. For physically 
interesting parameter choices the perturbations turned out to be stable. Depending on the 
scale, the results coincide with those of the ACDM model until redshifts of the order of 
z « 4 for k = 0.1 or z « 1.5 for k = 0.01. For z < 3 (k = 0.1) or for z < 0.7 (k = 0.01) the 
differences are substantial. On smaller scales there appear oscillations of the perturbation 
quantities similar to those in Chaplygin-gas cosmologies. For the moment we must leave open 
whether these oscillations represent a real effect inside the dark sector, potentially detectable 
by weak gravitational lensing and/or modifications of the CMB anisotropy spectrum, or 
whether they reflect a deficiency of the model, possibly due to the neglect of intrinsic non- 
adiabatic perturbations of the DE component. 

Acknowledgments 

The authors would like to thank Oliver Piattella for helpful discussions. This work 
was supported by the "Comision Nacional de Ciencias y Tecnologfa" (Chile) through the 
FONDECYT Grant No. 1110230 (SdC) and No. 1090613 (RH and SdC). J.F. and W.Z. 
acknowledge support by "FONDECYT-Concurso incentivo a la Cooperation Internacional" 
No. 1090613 as well as by CNPq (Brazil) and FAPES (Brazil). 



[1] A. G. Riess et al., Astron. J. 116, 1009 (1998); S. J. Perlmutter et al., Astrophys. J. 517, 565 
(1999); A. G. Riess et al., Astrophys. J. 607, 665 (2004); P. Astier et al., Astron. Astrophys. 
447, 31 (2006). 

[2] G. Hinshaw et al.. larXiv:1212.52261 

[3] S. Weinberg, Rev. Mod. Phys. 61, 1 (1989); E.J. Copeland, M. Sami, S. Tsujikawa, Int. J. 

Mod. Phys. D 15, 1753 (2006). 
[4] A. G. Cohen, D.B. Kaplan and A.E. Nelson, Phys. Rev. Lett. 82, 4971 (1999). 
[5] M. Li, Phys. Lett. B 603, 1 (2004). 
[6] S. D. H. Hsu, Phys. Lett. B 594, 13 (2004). 

22 



[7] D. Pavon and W. Zimdahl, Phys. Lett. B 628, 206 (2005). 

[8] W. Zimdahl and D. Pavon, Class. Quantum Grav. 24, 5641 (2007); W. Zimdahl, IJMPD 17, 
651 (2008). 

[9] Hsien-Chung Kao, Wo-Lung Lee and Feng-Li Lin, Phys. Rev. D 71, 123518 (2005); K. Enquist, 
S. Hannestad and M.S. Sloth, JCAP 0202, (2005) 004; Yungui Gong, Bin Wang and Yuan- 
Zhong Zhang, Phys. Rev. D 72, 043510 (2005); Xin Zhang and Feng-Quan Wu, Phys. Rev. 
D 72, 043524 (2005); Phys. Rev. D 76, 023502 (2007); Qiang Wu, Yungui Gong, Anzhong 
Wang andJ.S. Alcaniz, Phys. Lett. B 659, 34 (2008); Shao-Feng Wu, Peng-Ming Zhang, and 
Guo-Hong Yang, Class. Quantum Grav. 26 (2009) 055020; Yinzhe Ma, Yan Gong and Xuelei 
Chen, llarXiv:0901 .12151 

[10] R. Brustein and G. Veneziano, Phys. Rev. Lett. 84, 5695 (2000). 

[11] C. Gao, F. Q. Wu, X. Chen and Y. G. Shen, Phys. Rev. D 79, 043511 (2009). 

[12] R-G. Cai, B. Hu and Y. Zhang, Commun. Theor. Phys. 51, 954 (2009). 

[13] Lixin Xu and Yuting Wang, JCAP 1006, 002 (2010); iarXiv: 1006.02961 

[14] Xin Zhang, Phys. Rev. D79, 103509 (2009); arXiv:0901.22621 

[15] Rong-Jia Yang, Zong-Hong Zhu and Fengquan Wu, Int. J. Mod. Phys. A 26, 317 (2011); 
IarXiv: 1101 .47971 

[16] Chao-Jun Feng and Xin-Zhou Li, Phys. Lett. B 680, 184 (2009); l arXiv:0904.2"972l 

[17] K. Karwan and T. Thitapura, JCAP 1201, 017 f2012 ): la7XTv:1110.24511 

[18] Yuting Wang, Lixin Xu and Yuanxing Gui, Phys. Rev. D 84, 063513 (2011); arXiv:1110.440Tl 

[19] Yi Zhang and Hui Li, IarXiv: 1003.27881 

[20] I. Duran and D. Pavon, Phys. Rev. D 83, 023504 (2011); llarXiv:1012.2986l 
[21] Tian-Fu Fu, Jing-Fei Zhang, Jin-Qian Chen and Xin Zhang, Eur. Phys. J. C 72, 1932 (2012); 
IarXiv: 11 12.23501 

[22] Zhenhui Zhang, Song Li, Xiao-Dong Li, Xin Zhang and Miao Li, JCAP 1206, 009 (2012); 
IarXiv: 1204.61351 

[23] L.P. Chimento and M.G. Richarte, Phys. Rev. D 85, 127301 (2012); arXiv: 1207. 14921 
[24] L.P. Chimento, M. Forte and M.G. Richarte. larXiv:1301.2737l 
[25] B. Broda, IJMPD 21, 1250053 (2012): larXTv:1111.5785L 

[26] S. del Campo, J.C. Fabris, R. Herrera and W. Zimdahl, Phys. Rev. D 83, 123006 (2011). 
[27] J. Simon, L. Verde and R. Jimenez, Phys. Rev. D 71, 123001 (2005). 

23 



[28] D. Stern, R. Jimenez, L. Verde, M. Kamionkowski and S. A. Stanford, JCAP 1002, 008 
(2010). 

[29] R. Jimenez, L. Verde, T. Treu and D. Stern, Astrophys. J. 593, 622 (2003). 
[30] T-J. Zhang, C. Ma and T. Lan, Adv. in Astron. 2010, 184284 (2010). 
[31] C. Ma and T-J. Zhang, Astrophys. J. 730, 74 (2011). 

[32] O. Farooq, D. Mania and B. Ratra, Hubble parameter measurement constraints on dark energy, 
larXiv:1211.4253L 

[33] B.L. Lago, M.O. Calvao, S.E. Joras, R.R.R. Reis, I. Waga and R. Giostri, Type la su- 
pernova parameter estimation: a comparison of two approaches using current datasets, 
larXiv:1104.28741 

[34] R. Amanullah, et al. Astrophys. J. 716, 712( 2010). 

[35] D.L. Eisenstein et al. (SDSS) Astrophys. J. 633, 560 (2005). 

[36] C. Blake et all, The WiggleZ Dark Energy Survey: mapping the distance-redshift relation with 

baryon acoustic oscillations, arXiv:1108.2635 . 
[37] H.B. Sandvik, M. Tegmark, M. Zaldariaga and I. Waga, Phys. Rev. D 69, 123524 (2004). 
[38] W.S. Hipolito-Ricaldi, H.E.S. Velten and W. Zimdahl, Phys. Rev. D 82, 063507 (2010). 



24 



